Localized defect modes in graphene 
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We study the properties of localized vibrational modes associated with structural defects in a sheet of 
graphene. For the example of the Stone- Wales defects, one- and two-atom vacancies, many-atom linear va- 
cancies, and adatoms in a honeycomb lattice, we demonstrate that the local defect modes are characterized by 
stable oscillations with the frequencies lying outside the linear frequency bands of an ideal graphene. In the 
frequency spectral density of thermal oscillations, such localized defect modes lead to the additional peaks from 
the right side of the frequency band of the ideal sheet of graphene. Thus, the general structure of the frequency 
spectral density can provide a fingerprint of its quality and the type of quantity of the structural defect a graphene 
sheet may contain. 



I. INTRODUCTION 

Over the past 20 years nanotechnology has made an im- 
pressive impact on the development of many fields of physics, 
chemistry, medicine, and nanoscale engineering [1]. After the 
discovery of graphene as a novel material for nanotechnol- 
ogy [2], many properties of this two-dimensional object have 
been studied both theoretically and experimentally [3, 4]. 

Usually, depending on the type of growth and fabrication 
procedure employed, the graphene surface has structural de- 
fects which may change substantially its mechanical, elec- 
tronic, and transport properties. Structural defects may also 
change both electronic and phononic spectra of the graphene. 
Moreover, defects creates additional sources of scattering for 
phonons and electrons, so they may alter substantially the 
transport properties of graphene structures. 

By adding defects to an ideal graphene, we may change its 
properties, and this is the main concept of nanoengineering 
of defect structures on graphene. In particular, local defects 
may change the absorbing properties of graphene [5], whereas 
extended defects allow the use of graphene as a membrane 
for gas separation [6], and arrays of defects may change the 
conducting properties of graphene structures [7, 8]. 

We may divide various locaUzed defects in graphene into 
two classes: point defects (such as Stone-Wales defect, va- 
cancies, and adatoms) and one-dimensional dislocation-like 
defects [9]. Point defects in graphene act as scattering centers 
for electron and phonon waves [ 10-12] . Linear defects are re- 
sponsible for plastic deformations of graphene nanoribbons. 

In this paper, we study the oscillatory phonon modes lo- 
calized at the point defects in graphene. To the best of our 
knowledge, this problem was never analyzed in detail. In 
particular, we apply the molecular-dynamics numerical sim- 
ulations and demonstrate that each type of structural defects 
supports a number of localized modes with the frequencies 
located outside the frequency bands of the ideal graphene lat- 
tice. We also analyze nonlinear properties of these localized 
modes and demonstrate that their frequencies decrease when 
the input power grows. These localized defect modes manifest 
themselves in the thermalized dynamics of the graphene sheet, 
and they may be employed as fingerprints of the graphene 
structural quality. 



The paper is organized as follows. In Sec. II we describe 
our microscopic model and introduce the interaction poten- 
tials. Section III summarizes the types of the local defects in 
graphene analyzed in the paper. Sections IV and V are de- 
voted to the analysis of the properties of in-plane and out-of- 
plane defect modes. In Sec. VI and VII we analyze how the 
defect modes are manifested in the dynamics and thermalized 
relaxation of graphene structures. Section VIII concludes the 
paper. 



II. MODEL 

To model numerically the oscillations of the graphene sheet 
with local defects, we consider a rectangular graphene lattice 
of a finite extent (a graphene flake) with the size 10.8 x 10.6 
nm^ composed of TV = 4400 carbon atoms. We place a defect 
at the center of this graphene flake, and introduce periodic 
boundary condition to avoid the interaction with the bound- 
aries. 

To describe graphene oscillations, we write the system 
Hamiltonian in the form 
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where M ~ Ylniy is the mass of carbon atom {m^ is the 
proton mass), u„ = [x„(i),j/„(t), 2„(t)] is the radius vector 
of the carbon atom with the index n at the moment t. The 
term P„ describes the interaction of the atom with the index 
n and its neighboiing atoms. The potential energy depends 
on variations in bond length, bond angles, and dihedral angles 
between the planes formed by three neighboiing carbon atoms 
and it can be written in the form 

P = ^C/l+^C/2 + ^C/3 + ^C/4 + ^C/5, (2) 
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where 17^, with t = 1,2,3,4,5, are the sets of configura- 
tions including up to nearest-neighbor interactions. Owing 
to a large redundancy, the sets only need to contain configura- 
tions of the atoms shown in Fig. 1, including their rotated and 
mirrored versions. 
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FIG. 1 : Configurations containing up to ith nearest-neighbor interac- 
tions for (a) i = 1, (b) i = 2, (c) i = 3, (d) i = 4, and (e) i = 5. 



The potential [/i(uq,u^) describes the deformation en- 
ergy due to a direct interaction between pairs of atoms with 
the indexes a and (3, as shown in Fig. 1(a). The potential 
U2{ua, u^, u-y) describes the deformation energy of the an- 
gle between the valent bonds u^u^ and U/3U-y, see Fig. 1(b). 
Potentials [/^(uq,, u^, u-y, u^), i = 3, 4, 5, describes the de- 
formation energy associated with a change of the effective an- 
gle between the planes Uq,, u^, u-^ and u^, u^, u^, as shown 
in Figs. l(c,d,e)- 

We use the potentials employed in the modeling of the dy- 
namics of large polymer macromolecules [13-17] for the va- 
lent bond coupling, 

C/l(ui,U2) = ei{cxp[-Qo(p-/5o)] - 1}^, P= |U2 -uil, 

(3) 
where ei — 4.9632 eV is the energy of the valent bond and 
po = 1.418 A is the equilibrium length of the bond; the po- 
tential of the valent angle 

[/2(U1,U2,U3) = e2(cOS(/3-COS<po)^, (4) 
COSip= (U3 -U2,Ui -U2)/(|U3 - U2 1 ■ |U2 -Ui|), 

SO that the equilibrium value of the angle is defined as 
cosLpo — cos(27r/3) = ^1/2; the potential of the torsion 
angle 

C/j(ui,U2,U3,U4) = 6^(1 - ZiCOSCJ)), (5) 

008 0= (vi,V2)/(|vi| • |V2|), 
Vl = (U2 - Ui) X (U3 - U2), 
V2 = (U3 - U2) X (U3 - U4), 

where the sign z, = 1 for the indices i = 3, 4 (equilibrium 
value of the torsional angle (pa = 0) and Zi = —1 for the 
index i = 5 ((po = tt). 

The specific values of the parameters are ao ~ 1.7889 A^^, 
£2 — 1.3143 eV, and £3 = 0.499 eV, and they are found from 
the frequency spectrum of small-amplitude oscillations of a 
sheet of graphite [18]. According to the results of Ref. [19] 
the energy £4 is close to the energy £3, whereas £5 ^ £4 
(|£5/£4| < 1/20). Therefore, in what follows we use the val- 
ues £4 = £3 = 0.499 eV and assume £5 = 0, the latter means 
that we omit the last term in the sum (2). 

More detailed discussion and motivation of our choice 
of the interaction potentials (2), (3), (4) can be found in 
Ref. [20]. Such potentials have been employed for modeling 
of thermal conductivity of carbon nanotubes [21, 22] graphene 
nanoribbons [20] and also in the analysis of their oscillatory 
modes [23-25]. 



III. TYPES OF LOCALIZED DEFECTS IN GRAPHENE 

First, we introduce a model of graphene with defects. We 
take a rectangular flake of graphene with an ideal honeycomb 
lattice and, depending on the type of defect we wish to cre- 
ate, we add or remove some carbon atoms at its center by 
cutting the corresponding bonds. In this way, we define the 
atomic configuration {un}'^^-^ that will relax to the structure 
with a specific type of defect. To find the ground state of the 
graphene layer with defects, we should find the energy mini- 
mum for the interaction energy. 



N 

E 



P„ -> min : {u„}„^;^ 



(6) 



Problem (6) is solved numerically by means of the conju- 
gate gradient method. If {u'^}n=i is the ground state of this 
rectangular flake of graphene, then for small-amplitude oscil- 
lations we can write u„(i) = ul\ + v„(f), where |v„| ^ po- 
Then, the equations of motion corresponding to the Hamilto- 
nian (1) can be written as a system of 3A^ linear equation for 
3iV variables. 
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(7) 
To find all linear modes of the graphene sheet, we need to 
find numerically all 3iV x 3A^ eigenvalues and corresponding 
eigenvectors of the real symmetric matrix B = {Bjn)^n=i- 
If we define A and e = {v"}^^]^ as the eigenvalue and 
normalized eigenvector, namely Be = Ae and (e, e) — 
^jj(v",v") = 1, then the solution of Eq. (7) will have the 
form v„(i) = Av° exp(iwt), where lj = ^J\/M) is the 
mode frequency and A is the mode amplitude. 

We characterize the degree of spatial localization of any 
of the oscillatory eigenmode by the parameter of localization 



(inverse participation number), d = ^ 



^ --",v0)2. For 
the modes which are not localized in space, d « 1/A^, for 
the mode localized on a single atom, d = \. Inverse value 
Nd = l/d (participation number) characterizes the number of 
atoms which are involved into this oscillatory mode. 

For an ideal flake of graphene, all eigenmodes are delocal- 
ized and for the periodic boundary conditions they are dis- 
tributed homogeneously on the rectangular surface. All such 
oscillatory modes can be divided into two classes; (i) in- 
plane modes, when the atoms oscillate along the flat surface 
of the graphene sheet, and (ii) out-of -plane modes, when the 
atoms oscillate perpendicular to the surface of the graphene 
sheet. The frequencies of the former modes are located in 
the domain < w < 1600 cm~^, whereas the frequencies 
of the latter modes correspond to a narrow spectral range, 
<a; <900cm-i. 

Bellow, first we study the in-plane localized defects shown 
in Fig. 2 which do not bend the surface of the graphene sheet 
(Sec. IV), but then consider the case of out-of-plane defects 
(Sec. V). 




FIG. 2: Localized defects which do not change the flat surface of 
the graphene layer, (a) Stone-Wales defect SW(55-77), (b) single va- 
cancy Vi(5-9); (c) double vacancy V2(5-8-5); (d) double vacancy 
V2(555-777); (e) double vacancy V2(5555-6-7777); (f) quadruple 
vacancy V4(555-9); (g) vacancy Vg(5-7-66-7-5), extended linear de- 
fect obtained by removing a zigzag-like array of 8 carbon atoms; (h) 
vacancy Ve, a defect created by removing 6 carbon atoms with the 
formation of a hexagonal hole. Grey color is employed to show a 
change of the lattice structure introduced by the presence of a defect. 



IV. OSCILLATIONS OF IN-PLANE DEFECTS 

To create the Stone-Wales defect SW(55-77), we should ro- 
tate one valent bond by 90 degrees. Then four old neighbor- 
ing valent bonds break creating four new bonds. As a result, 
we create two pentagons and two heptagons (55-77) instead 
of four hexagons of the perfect honeycomb lattice, as shown 




FIG. 3: Energy of in-plane oscillations for the Stone-Wales defect 
SW(55-77) with the frequency ujd ~ 1613.2 cm~^, and the partici- 
pation number Nd ~ 26.31. 



in Fig. 2(a). In the ground state, this defect has the energy 
Ed = 2.06 eV, and it is stable. Although the defect's energy 
is higher than that of the ground state by the value Ed, the 
reverse change of the bonds requires to overtake the energy 
banierof A£' = 9eV[9]. 

Our analysis of the localized modes of the lattice with the 
defect shows that the Stone-Wales defect supports four in- 
plane localized modes with the frequencies tud = 1613.2, 
1614.1, 1615.6, 1616.3 cm^^ (and the participation numbers 
Nd = 26.31, 27.06, 25.65, 26.14) and also four out-of-plane 
localized modes with the frequencies ujd = 904.7, 907.5, 
909.4, 935.1 cm-i (Nd = 44.74, 21.81, 37.24, 7.90). The 
characteristic profile of the oscillatory energy of this defect 
mode is shown in Fig. 3, where we observe that the energy is 
localized entirely in the region of defect, so the used periodic 
boundary conditions do not have any influence on the mode's 
frequency and shape. The oscillatory dynamics of all local- 
ized modes is shown in movies from Supplementary Material. 

In order to create a single vacancy Vi(5-9) in the graphene 
lattice, we remove just one carbon atom breaking three va- 
lent bonds. The remaining lattice relax with two neighbor- 
ing atoms linked by a new valent bond, and instead of three 
hexagons the structure (5-9) is formed with one pentagon and 
one nonagon, see Fig. 2(b). In this defect structure, one car- 
bon atom has only two valent bonds. In order to keep the 
interaction potentials those bonds unchanged in our calcu- 
lations, we assume that this atom is coupled to a hydrogen 
atom and has mass M = 13m.p. Them the energy of the 
defect Vi(5-9) is found to be Ed = 2.26 eV. This defect sup- 
ports three in-plane localized eigenmodes with the frequen- 
cies uJd = 1600.6, 1601.3, 1606.3 cm^^ (Nd = 146.99, 
79.41, 14.58) and two out-of-plane lightly localized modes 
with the frequencies cod = 1227.0, 1791.0 cm^^ (Nd = 5.03, 
2.65). 
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FIG. 4: Deformation of the graphene sheet in the region of location of 
(a) inverse Stone-Wales defect l2(7557) (two carbon adatoms with a 
valent bond are inserted in a perfect graphene lattice; they are shown 
by larger balls) and (b) two adatoms defect A2(5556777). Grey color 
is employed to show the change of the lattice structure introduced by 
the presence of a defect. 



Double-vacancy defect V2(5-8-5) shown in Fig. 2(c) can be 
created by removing two neighboring carbon atoms, break- 
ing four valent bonds. After this procedure, the remaining 
four neighboring carbon atoms create two new valent bonds, 
and in the lattice structure four hexagons are replaced by 
two pentagons and one octagon. The energy of the double- 
vacancy defect V2(5-8-5) is found to be Ed = 3.44 eV. This 
defect supports four in-plane eigenmodes of localized oscil- 
lations with the frequencies Ud = 1605.9, 1606.0, 1619.2, 
1619.3 cm-i (Nd = 65.91, 66. 16, 15.21, 15.32) and four out- 
of-plane localized modes with the frequencies ujd ~ 902.6, 
905.0, 986.4, 1015.0 cm^^ {Nd = 64.16, 57.77, 6.18, 7.27). 

In order to create the double-vacancy defect V2(555-777), 
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FIG. 5: Dependencies of the energy E of localized defect modes 
vs. frequency uj, for two types of defects: inverse Stone-Wales de- 
fect l2(7557), shown by blue lines, and for the Stone- Wales defect 
SW(55-77), shown by red lines. Parts (a) and (b) show out-of-plane 
and in-plane oscillations, respectively. Vertical dotted line marks the 
edge of the linear spectrum band of phonons of the ideal graphene 
sheet. 
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FIG. 6: Frequency spectral density of thermal oscillations of a 
graphene sheet with inverse Stone-Wales defect I2, when two 
adatoms are inserted into the honeycomb lattice. The graphene sam- 
ple consists of TV = 578 atoms and it has the size of 3.93 x 3.83 
nm^, being placed at temperature of T = 300K. Vertical solid bars 
show the frequencies of the localized oscillation of the defect modes. 



we should start from the defect V2 (5-8-5) and rotate one of 
the bonds of its octagon by 90 degrees. Then four old bonds 
break and four new bonds appear In the graphene lattice we 
obtain instead of seven hexagons three pentagons and three 
heptagons (555-777), as shown in Fig. 2(d). This defect is 
characterized by the energy Ed = 2.91 eV. Our analysis of 
the eigenmodes localized on this defect demonstrate that it has 
no in-plane localized modes at all, but it supports four out-of- 
plane modes with the frequencies cud = 909.1, 916.0, 916.0, 
916.6 cm-i (iVd = 46.22, 15.87,15.87, 11.30). 

Double-vacancy defect V2(5555-6-7777) is formed from 
the defect V2(555-777) by another bond rotation. In the hon- 
eycomb lattice, instead of ten perfect hexagons we have four 
pentagons, one hexagons, and four heptagons, as shown in 
Fig. 2(e). The energy of this double-vacancy defect is Ed = 
2.96 eV. This defect supports seven localized in-plane eigen- 
modes with the frequencies ujd ~ 1605.5, 1605.6, 1607.3, 
1607.4, 1666.4, 1773.7, 2009.4 cm^^ (and the participation 
numbers Nd = 25.84, 25.50, 25.99, 23.65, 5.86, 4.57, 3.05) 



and four out-of-plane localized modes with the frequencies 
ojd = 905.7, 907.2, 916.4, 918.4 cm-^ {N^ = 46.72, 54.22, 
15.07, 14.52). 

The largest hole appears when we remove 4 neighboring 
carbon atoms from the graphene lattice with Y-shape va- 
lent bonds. After that, six carbon atoms coupled to the re- 
moved atom got involved into new valent bonds, thus creat- 
ing quadruple vacancy defect V4(555-9). In this defect, in- 
stead of six hexagons we have three heptagons and one large 
nonagon, as shown in Fig. 2(f). The energy of this defect is 
Ed = 5.37 eV. This defect supports three in-plane eigenmodes 
with the frequencies cjrf = 1601.9, 1602.4, 1602.4 cm^^ (par- 
ticipation numbers N^ = 59.02, 40.46, 40.21) and also four 
out-of-plane modes with the frequencies utd ~ 902.3, 908.9, 
976.6, 997.3 cm-i (Na = 52.27, 11.01,7.23,6.51). 

For comparison, next we consider an extended defect in the 
form of the vacancy V8(5-7-66-7-5), that appears after remov- 
ing a zigzag chain of an array of eight carbon atoms shown 
in Fig. 2(g). In the honeycomb lattice, instead of ten per- 
fect hexagons we create two pentagons, two hexagons, and 
two heptagons; the defect energy is Ed ~ 11.92 eV. This ex- 
tended defect supports seven in-plane localized modes with 
the frequencies Wd = 1604.2, 1606.8, 1607.2, 1613.4, 1624.3, 
1627.3, 1653.2 cm"! (Nd = 57.72, 19.39, 18.03, 11.53, 
10.67, 1 1 .49, 4.36) and also four out-of-plane localized modes 
with the frequencies ujd = 909.3, 910.5, 927.2, 1091.0 cm^^ 
corresponding to Nd — 21.35, 21.53, 6.25, and 2.69. 

If we remove one hexagon from the honeycomb lattice, we 
create a vacancy Vg that can not be "healed" by the lat- 
tice transformations. This vacancy has the form of a perfect 
hexagon hole with the diameter of 0.57 nm, with almost un- 
changed old valent bonds, see Fig. 2(h). Our analysis demon- 
strates that such holes do not support any localized modes. 



V. OSCILLATIONS OF OUT-OF-PLANE DEFECTS 

Unlike vacancies, adding extra atoms to the graphene sheet 
lead to a deformation of the planar structure of the lattice. In 
the place where we have more atoms, the graphene sheet is 
no longer flat, and it becomes either concave or convex, trans- 
forming into one of the two equivalent non-planar states. 

When two migrating adatoms meet each other and form a 
dimer, they can be incorporated into the graphene lattice, cre- 
ating the so-called inverse Stone-Wales defect l2(7557). In 
such a state, instead of three hexagons we have two pentagons 
and two heptagons, as shown in Fig. 4(a). The surface of the 
sheet become deformed bending into either side. 

Our analysis shows that such a defect supports 8 localized 
modes for which atoms move along the curved surface. The 
mode eigenfrequencies are: ujd = 1604.8, 1613.1, 1624.2, 
1641.3, 1644.2, 1650.7, 1723.5, 1749.4 cm^^ (with the partic- 
ipation numbers Nd = 73.63, 5.93, 19.82, 6.73, 8.64, 12.37, 
3.14, 2.85). This defect also supports two out-of-plane local- 
ized modes for which atoms move perpendicular to the sur- 
face, the modes have the frequencies ojd = 902.4, 903.9 cm^^ 
(A^d = 99.55, 71.79). 

An excess of two carbon atoms may lead to another concave 



deformation of the graphene lattice and two adatom defect 
mode A2 (5556777), when instead of 6 perfect hexagons we 
have 3 pentagons, 1 hexagon, and 3 heptagons, see Fig. 4(b). 
Such a defect may appear when we remove 4 carbon atoms 
from the graphene lattice with Y-shape valent bonds but insert 
instead of them 6 carbon atoms with a structure of a perfect 
hexagon. The transverse oscillations of this defect blister is 
considered in Ref. [26]. As can be seen from Fig. 4, the defect 
I2 has a higher density of the atom localization in comparison 
with the defect A2, and as a result the defect I2 will bend the 
graphene sheet stronger than the defect A2. 

Our analysis of the localized oscillations of the graphene 
sheet demonstrate that the defect A2 (5556777) supports 11 
localized modes for which the atoms move along a concave 
surface of the graphene sheet. The frequencies of those os- 
cillations are: ujd = 1601.0, 1603.3, 1604.2, 1608.6, 1610.7, 
1613.6, 1614.6, 1617.8, 1682.4, 1704.2, 1706.4 cm"! (par- 
ticipation numbers Nd = 131.01, 19.52, 22.49, 13.37, 7.28, 
7.96, 19.94, 4.48, 3.94, 3.52, 2.77). This defect supports also 
two out-of-plane oscillatory modes for which atoms move per- 
pendicular to the graphene sheet; these modes have the fre- 
quencies: ojd = 900.7, 910.4 cm-i {Nd = 85.54, 6.15). 



VI. MODELING OF LOCALIZED OSCILLATIONS 

To verify the results of our modal analysis, we model the 
oscillatory dynamics of a graphene with defects. We con- 
sider a graphene sheet of a finite extent, with the size of 
10.8 X 10.6 nm^ with a local defect at its center To model the 
energy spreading in an effectively infinite graphene sheet, we 
introduce lossy boundary conditions assuming that all edge 
atoms experience damping with relaxation time t,. = 10 ps. If 
the excitation of the localized eigenmode lead to the creation 
of a localized state, or breathing-like mode, the kinetic energy 
not vanish but instead it will approach a certain nonzero value. 
Otherwise, the initial excitation will vanish completely, so that 
the kinetic energy 



Ek 
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A/(u„,u„) 



will vanish when t — !> 00. 

Our numerical simulation results demonstrate that all lo- 
calized modes with the frequencies located outside the lin- 
ear frequency bands of the ideal graphene sheet are stable. 
For small amplitudes, when the energy of the eigenmode less 
than 0.06 eV, the amplitude and frequency of the mode do not 
change, and this result indicates that in that limit the oscilla- 
tions are described by a linear theory. However, when more 
energy is pumped into the mode, we observe nonUnear effects 
when the frequency decreases with the energy, as shown in 
Fig. 5. In that case, the energy of the defect mode may exceed 
leV. 

We notice that the diagonalization of the matrix of the 
second derivatives B = {Bjn)fn=i ^^^^ ^^^ &^^ localized 
modes with the frequencies within the frequency bands of 
the ideal graphene sheet. However, numerical simulations 



demonstrate that such modes are not stable and the amplitude 
decays very rapidly. 



VII. INFLUENCE OF DEFECTS ON THE GRAPHENE'S 
FREQUENCY SPECTRA 

Our results summarized above suggest that all frequencies 
of localized defect modes are located outside the linear bands 
of the oscillation spectrum of an ideal graphene sheet. This 
result is expected, because any oscillation from the continu- 
ous spectrum is not spatially localized, and it corresponds to 
one of the extended modes. The splitting of the continuous 
modes of a perfect sheet of graphene, and localized modes 
of the graphene with defects can be employed as an impor- 
tant fingerprint of the quality of the graphene honeycomb lat- 
tice. We demonstrate this approach for the case of a graphene 
sheet of a finite extent of 3.93 x 3.83 nm^ with inverse Stone- 
Wales defect I2. This patch of graphene has 578 carbon atoms 
two of which create one point-like defect of the type shown in 
Fig. 4(a). 

In order to find the frequency spectrum of thermal oscilla- 
tions of a graphene sheet, we place it in a Langevin thermostat 
at temperature T = 300 K. After complete thermalization, we 
disconnect the graphene from the thermostat and study the os- 
cillation spectra of the excited thermal oscillations. 

Frequency spectral density of thermal oscillations of a sheet 
of graphene is shown in Fig. 6. As follows from those re- 
sults, several additional peaks appear in the frequency spec- 
trum above the upper edge of the spectral band of linear oscil- 
lations of a perfect graphene sheet (w > 1600 cm^^). These 



peaks correspond to the eigenmodes of the localized oscilla- 
tions of the defect I2. Thus, from the specific structure of the 
frequency spectral density we may judge about the types and 
density of structural defects in a particular graphene sample. 



VIII. CONCLUSIONS 

We have studied the effect of localized structural defects on 
the vibrational spectra of the graphene's oscillations. We have 
demonstrated that the local defects (such as Stone- Wales de- 
fects, one- and two-atom vacancies, many-atom linear vacan- 
cies and adatoms) are characterized by stable localized oscil- 
lations with the frequencies lying outside the linear frequency 
bands of the ideal graphene. In the frequency spectral den- 
sity of thermal oscillations, such localized defect modes lead 
to the additional peaks located on the right side from the fre- 
quency band of the ideal sheet of graphene. Thus, the general 
structure of the frequency spectral density can provide a fin- 
gerprint of its quality and the type of quantity of the structural 
defect a graphene sheet may contain. 



Acknowledgements 

Alex Savin acknowledges a warm hospitality of the Non- 
linear Physics Center at the Australian National University, 
and thanks the Joint Supercomputer Center of the Russian 
Academy of Sciences for the use of their computer facilities. 
The work was supported by the Austrahan Research Council. 



[1] K.E. Drexler, Nanosystems: Molecular Machinery, Manufac- 
turing, and Computation (New York, Wiley, 1992). 

[2] K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, Y Zhang, 
S.V. Dubonos, I.V. Grigorieva, and A. A. Firsov, Science 306, 
666 (2004). 

[3] A.K. Geim and A.H. MacDonald, Phys. Today 60, 35 (2007). 

[4] A.H. Castro Neto, F. Guinea, N.M. Peres, K.S. Novoselov, and 
A.K. Geim, Rev. Mod. Phys. 81, 109 (2009). 

[5] Y.-H. Zhang, K.-G. Zhou, K.-F Xie, X.-C. Gou, J. Zeng, H.- 
L. Zhang, and Y Peng, J. Nanosci. Nanotechnology, 10, 7347 
(2010). 

[6] X. Qin, Q. Meng, Y Feng, and Y Gao, Surface Science 607, 
153 (2013). 

[7] D.J. Appelhans, L.D. Carr and M.T. Lusk, New J. Phys. 12, 
125006 (2010). 

[8] J. Song, H. Liu, H. Jiang, Q. Sun, and X. C. Xie, Phys. Rev B 
86, 085437 (2012). 

[9] F. Banhart, J. Kotakoski, and A. V. Krashenirmikov, ACS Nano 
5,26(2011). 
[10] N. Gorjizadeh, A. A. Farajian and Y Kawazoe, Nanotechnol- 
ogy 20, 015201 (2009). 
[11] J.-H. Chen, W. G. Cullen, C. Jang, M. S. Fuhrer, and E. D. 

Williams, Phys. Rev. Lett. 102, 236805 (2009). 
[12] J. Haskins, A. Kinaci, C. Sevik, H. Sevincli, G. Cuniberti, and 

T. Cagin ACS Nano 5, 3779 (201 1). 
[13] D. W. Noid, B. G. Sumpter, and B. Wunderlich, Macro- 



molecules 24, 4148 (1991). 
[14] B. G. Sumpter, D. W. Noid, G. L. Liang, and B. Wunderlich, 

Adv. Polym. Sci. 116, 27 (1994). 
[15] A.V. Savin and L.I. Manevitch, Phys. Rev. B 58, 11386 (1998). 
[16] A.V. Savin and L.I. Manevitch, Phys. Rev. E 61, 7065 (2000). 
[17] A. V. Savin and L. I. Manevitch, Phys. Rev. B 67, 144302 

(2003). 
[18] A. V. Savin and Yu. S. Kivshar, Europhys. Letters 82, 66002 

(2008). 
[19] D. Gunlycke, H. M. Lawler, and C. T. White, Phys. Rev. B 77, 

014303 (2008). 
[20] A. V. Savin, Yu. S. Kivshar, and B. Hu. Phys. Rev. B 82, 195422 

(2010). 
[21] A. V. Savin, Yu. S. Kivshar, and B. Hu. Europhys. Lett. 88, 

26004 (2009). 
[22] A.V. Savin, B. Hu, and Yu. S. Kivshar. Phys. Rev. B 80, 195423 

(2009). 
[23] A. V. Savin and Yu. S. Kivshar. Appl. Phys. Lett. 94, 111903 

(2009). 
[24] A. V. Savin and Yu. S. Kivshar. Europhys. Lett. 89, 46001 

(2010). 
[25] A. V. Savin and Yu. S. Kivshar. Phys. Rev. B 81, 165418 (2010). 
[26] C.-W. Pao, T.-H. Liu, C.-C. Chang, and D.J. Srolovitz, Carbon 

50, 2870 (2012). 



This figure "I2_902.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_903.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_904.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_905.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_905.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_906.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_907.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_908.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555777_909.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_909.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_910.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_910.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_916.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555777_916a.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555777_916b.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555777_917.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_918.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_927.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_935.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_986.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_1015.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1089.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "Vl_1227.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "Vl_1600.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "Vl_1601.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1605.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1605a.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1605b.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1605.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "Vl_1606.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_1606.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1607a.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1607b.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1607.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1608.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1613.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_1613.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_1614.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1614.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_1615.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "SW5577_1616.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2585_1619.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1624.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1624.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1628.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1641.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1644.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1651.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V8_1654.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1666.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1724.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "I2_1749.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_1773.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "Vl_1791.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



This figure "V2555567777_2009.gif" is available in "gif ' format from: 



http://arxiv.org/ps/1304.7702vl 



